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I.   Introduction 

The  family  of  models  of  computation  used  In  this  paper  Is  the  parallel 
random-access-machines  (PRAMs).  All  members  of  this  family  employ  p  synchronous 
processors  all  having  access  to  a  common  memory.  We  mention  3  members  of  the  PRAM 
family  In  descending  order  of  strength.  In  a  concurrent-read  concurrent-write 
(CRCW)  PRAM  simultaneous  reading  from  the  same  memory  location  is  allowed  as  well 
as  simultaneous  writing.  In  the  latter  case  the  lowest  numbered  processor 
succeeds.  A  concurrent-read  exclusive-write  (CREW)  PRAM  allows  simultaneous 
reading  into  the  same  memory  location  but  not  simultaneous  writing.  An  EREW  PRAM 
does  not  allow  simultaneous  reading  or  writing.  See  [Vl-83]  for  a  recent  survey 
of  results  concerning  the  PRAM  family.     .    -    . 

Recall  that  a  bridge  in  a  connected  graph  Is  an  edge  whose  removal 
disconnects  the  graph  and  a  digraph  is  strongly  connected  if  for  every  two 
vertices  u  and  v  there  is  a  directed  path  from  u  to  v.  The  following  problem  is 
considered.  Input.  A  connected  bridgeless  undirected  graph  G-(V,E),  |v|-n, 
|E|=m.  The  strong  ortentatton  problem.  Assign  directions  to  its  edges  so  that 
the  resulting  digraph  is  strongly  connected.  It  is  known  that  such  strong 
orientation  exists  iff  the  input  graph  is  connected  and  bridgeless.  [At-8A] 
indicates  that  this  problem  has  a  simple  linear  time  serial  solution:  Apply 
depth-first  search  to  the  graph.  Orient  tree  edges  away  from  the  root  and  back 
edges  towards  the  root.  Since  depth-first  search  seems  to  be  Inherently  serial 
(see  (R-83]  for  some  evidence)  another  approach  is  required  for  designing 
efficient  parallel  algorithms.  [At-BA]  gave  an  O(log  n)  time  algorithm  for  this 
problem  using  n^  processors  on  a  CRCW  PRAM.  Atallah's  algorithm  runs  in 
0(n^/p  +  log^n)  time  using  p  processors  on  a  CREW  PRAM.  Our  algorithm  runs  in 
O(log  n)  time  using  only  n-Hn  processors  on  a  CRCW  PRAM.  An  alternative 
implementation  of  the  algorithm  runs  In  0(n^/p)  time  using  p  <  n  /log  n  processors 
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on  an  CREW  PRAM.  This  is  optimal  for  dense  graphs.  A  very  high-level  description 
of  our  algorithm  is  close  to  Atallah's  algorithm.  However,  only  a  complete  and 
drastic  change  in  the  way  parallelism  is  approached  allows  the  claimed 
improvements  in  efficiency.  Algorithmic  methods  that  reduce  the  number  of 
processors  without  changing  the  running  time  are  of  great  importance  for  the 
design  of  parallel  algorithms.  The  survey  [Vi-83]  mentions  many  parallel 
algorithm  that  achieved  a  previously  known  running  time  with  less  processors. 
[TV-83]  introduced  a  useful  technique  for  computing  various  functions  on  trees 
using  an  Euler  path  that  visits  each  edge  of  the  tree  exactly  twice.  Substantial 
simplifications  and  amlif Ications  of  this  technique  are  also  given.  This 
technique  needs  O(log  n)  time  using  0(n)  space  and  n  processors  on  an  EREW  PRAM. 

We  assume  throughout  that  our  input  graph  is  connected  and  bridgeless.  If  we 
do  not  know  it  we  can  check  connectivity  using  the  algorithms  of  [SV-82]  or 
[CLC-81]  or  [Vi-81]  depending  on  the  complexity  efficiency  being  sought  (O(log  n) 
time  using  n+m  processors  or  the  alternative  one).  Existence  of  bridges  can  be 
checked  using  the  algorithm  of  [TV-83]  within  the  efficiencies  claimed  here. 

In  order  to  keep  this  presentation  short  our  exposition  will  focus  on  the 
main  algorithmic  aspects  while,  like  [At-841,  tedious  implementation  details  will 
be  suppressed.  Also,  there  are  several  cases  where  a  more  involved  implentatlon 
would  enable  big  portions  of  the  algorithm  to  run  on  an  EREW  PRAM  instead  of  the 
CREW  PRAM  as  presented. 

IIv  A  very  high-level  desGrtption  of  the  algorithm 

Let  T  be  an  undirected  spanning  tree  of  G  and  H  be  a  directed  version  of  T 
rooted  at  vertex  r. 

(a)  Find  T  and  H.  Let  q=lE-T|.  Assign  serial  numbers  from  1  to  q  to  edges  of 
E-T,  so  that  E-T  =  {e.,...,e  }  .   Remark.    [At-8A]   allowed   this   assignment   of 
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serial  numbers  to  edges  of  E-T  to  be  done  In  any  order.  However,  our  algorithm 
introduces  a  refinement  into  this  assignment  which,  perhaps  surprisingly,  turns 
out  to  be  very  powerful  since  it  enables  the  improvements  in  efficiency. 
Let  C^  be  the  fundamental  cycle  which  is  Induced  by  ej^,  1  <  k  <  q  .  Let  a  be  an 
edge  of  T  and  C,  Ci  ,...,Cj  be  all  fundamental  cycles  that  contain  a,  where 
1  <  i^  <  ±2  <  ...  <  i-i  <  (\  -  Define  the  lowest  numbered  cycle,  C^  ,  to  be  the 
master  cycle  of  a  and  e^^   to  be  the  master  edge  of  a  (denoted  MASTER(a)). 

(b)  Assign  an  arbitrary  direction  to  every  edge  in  E-T.  To  each  a  in  T 
assign  the  direction  which  conforms  with  the  direction  of  MASTER(a)  in  the  master 
cycle  of  a .        " 

Theorem  ( { A t-8 A ] ) ♦  The  algorithm  produces  a  strong  orientation  of  G. 

III.-  The  algorithm 

The  algorithm  presented  below  has  two  implementations.  One  runs  in  O(log  n) 
time  using  n+m  processors  on  a  CRCW  PRAM.  The  other  runs  in  0(n2/p)  time  using  p 
processors  on  a  CREW  PRAM.  -  "'       '  '  - 

Step  -1.  A  spanning  tree  T  of  G  is  computed.  For  this  we  use  an  adaptation 
of  the  connectivity  algorithm  of  [SV-82]  as  described  in  [AS-83]  and  [TV-83].  It 
runs  in  O(log  n)  time  using  n+m  processrs  on  a  CRCW  PRAM.  The  alternative 
Implementation  uses  the  0(n2/p)  time  algorithm  of  [CLC-81]  that  uses  p  <  n^/log^n 
processors  on  a  CREW  PRAM. 

The  efficiency  of  the  rest  of  the  algorithm  is  optimal  for  graphs  with 
m  =  f2(nlog  n):  It  takes  O(log  n)  time  using  (nlog  n  +  m)/log  n  processors  on  a 
CREW  PRAM.  So  if  it  will  be  possible  in  the  future  to  find  a  spanning  tree  within 
this  optimal  efficiency  only  Step  1  will  have  to  change. 

Step  -2.  Compute  H,  a  directed  version  of  T  which  is  rooted  at  some  vertex  r. 
We  summarize  below  a  solution  of  [TV-83]  for  step  2  and  3  since  extentions  of  this 
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solutlon  are  needed  for  later  steps  of  the  present  algorithm.  (a)  Replace  each 
edge  (u,v)  of  T  by  two  anti-parallel  directed  edges  u  ->^  v  and  v  -»^  u.  The  new 
digraph  is  called  T'.  In  (d)  below  we  determine  for  each  edge  e  in  T'  whether  it 
(or  its  anti-parallel  edge)  is  oriented  away  from  the  root  r,  thereby  determining 
H.  Since  in-degree(v)  =  out-degree(v)  for  every  vertex  v  in  T',  T'  has  an  Euler 
path  that  starts  and  ends  at  the  root  r.  Next,  we  compute  this  path  into  a  vector 
named  D.  (b)  For  each  vertex  v  of  T  we  do  the  following.  Let  degree(v)  =  d  in  T 
and  let  the  d  adjacent  edges  of  v  in  T  be  (v,Uj ),...,  (v.u^j) .  Set 
D(uj.  -^  v)  :=  V  -^  Uj^j  jjjqJ  j  for  1  <  1  <  d.  Now  D  has  an  Euler  cycle.  The 
"correction"  D(uj  ■*■  r)  :=  'end  of  list'  (where  d  =  degree(r))  gives  the  Euler 
path.  Next,  the  distance  of  each  edge  of  T'  from  the  end  of  the  path  is  computed 
Into  a  vector  R  using  a  'doubling'  procedure,  (c)  Initialize:  R(e)  :=  1  for  all 
edges  e  in  T',  and  R('end  of  list')  ":=  0.  Apply  flog  2n-2l  iterations  in  parallel 
(2n-2  is  the  length  of  the  Euler  path):  R(e)  :=  R(e)  +  R(D(e));  D(e)  :-  D(D(e)). 
(d)  if_  R(e)  >  R( 'anti-parallel  edge  of  e')  then  e  (and  not  its  anti-parallel  edge) 
is  oriented  away  from  the  root  r.  Coraplextty.  The  computation  of  (c)  (resp.  (a), 
(b)  and  (d))  needs  O(log  n)  (resp.  0(1))  time  using  n  processors.  Remarks.  (1) 
Let  f  be  an  edge  in  T.  For  better  understanding  of  some  of  the  following  steps  of 
the  algorithm  try  to  think  about  the  (directed)  copy  of  f  e  H  in  the  Euler  path  as 
a  left  parenthesis  and  its  anti-parallel  edge  as  a  right  parenthesis.  If  we  do  so 
for  every  such  f  the  Euler  path  will  correspond  to  a  legal  sequence  of 
parentheses,  where  matching  pairs  of  parentheses  will  represent  the  two  copies  of 
an  edge  of  T.  (2)  In  order  to  obtain  initially  a  "desired"  representation  of  the 
tree  [TV-83]  proposes  to  consider  application  of  the  sorting  algorithm  of  [AKS-83] 
to  the  edges  of  the  tree.  This  algorithm  runs  in  O(log  n)  time  using  n  processors 
on  an  EREW  PRAM.   We  avoid  elaborating  on  this. 

Step  3.    Compute   preorder  and  postorder  numbering  of  the  vertices  of  H  and 


their  levels  (length  of  their  directed  path  from  r  In  H).  For  all  edges  ec H 
(resp.  ecT'-H)  initialize  R(e)  :=  1  (resp.  R(e)  :=  0  ).  The  doubling  procedure 
assigns  to  R(e)  of  each  edge  eeT'  how  many  vertices  were  entered  for  the  first 
time  by  the  Euler  path  following  the  appearance  of  e  on  this  path.  It  is  a  simple 
exercise  to  compute  preorder  numbering  from  this  information.  Computation  of 
postorder  numbering  is  obtained  similarly:  For  all  edges  eeH  (resp.  eeT'-H) 
initialize  R(e)  :=  0  (resp.  R(e)  :=  1  ).  The  doubling  procedure  assigns  to  R(e) 
of  each  edge  eeT'-H  how  many  vertices  were  exited  for  the  last  time  by  the  Euler 
path  following  the  appearance  of  e  on  this  path.  And  apply  the  previous  exercise. 
Computation  of  levels.  For  all  edges  eeH  (resp.  eeT'-H)  initialize  R(e)  :=  1 
(resp.  R(e)  :=  -1  ).  Following  the  doubling  procedure  -R(e)  of  each  edge 
e  e  T'-H  is  the  level  of  both  e  and  its  head  vertex  in  H.  Complexity.  The 
computation  of  preorder,  postorder  and  the  levels  need  O(log  n)  time  using  n 
processors. 

Let  Uj,...,u£  be  all  vertices  which  are  adjacent  to  it  by  an  edge  of  E-T. 
Denote  by  Ad-LCA(v)  the  lowest  common  ancestor  in  H  of  v,Uj,...,U£  . 
(Ad-LCA(v)  =  LCA(v,Uj, ... . ,U£)  J.  We  say  that  an  edge  e=(v,u)  e  E-T  causes^ 
Ad-LCA(v)  if  LCA(v,u)  =  Ad-LCA(v). 

Step  4.1.  For  every  vertex  v  compute  Ad-LCA(v)  and  if  Ad-LCA(v)  is  different 
than  V  select  (arbitrarily)  an  edge  e  that  causes  it.  This  is  done  as  follows. 
Let  the  predicate  ANCESTOR(u,v)  be  "true"  iff  u  is  ancestor  of  v  in  H.  Recall  (see 
[Ta-7A])  that, 

ANCESTOR(u,v)  =  [PREORDER(u)  <  PREORDER(v)]AND[POSTORDER(u)  >  POSTORDER(v) ] . 
Partition  the  £  edges  (v,Ui ),.,., (v,uo )  of  E-T  which  are  adjacent  to  v   into   four 
sets.   Edges  where: 

(a)  [ P REORDER (uj^)  >  PREORDER(v)  ]  AND[POSTORDER(u^)  <  POSTORDER(v)  ]  .  In  this  case  v 
is  an  ancestor  of  u,  .   Therefore,  these  edges  cannot  affect  Ad-LCA(v). 
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(b)  [PREORDER(uj^)  <  PREORDER(v)]ANDIPOSTORDER(u^)  <  POSTORDER(v) ] .  Find  among 
these  edges  the  one  which  has  minimum  PREORDER(u,  ) .   Denote  this  edge  by  (v,uo  ). 

(c)  [PREORDER(u^)  >  PREORDER(v)]AND[POSTORDER(u^^)  >  POSTORDER(v)].  (1)  Find  among 
these  edges  the  one  which  has  maximum  POSTORDER(uj^) .  Denote  this  edge  by  (v.u^  ). 
(2)  Find  LCA(v,U£  )  and  LCA(v,U£  )  in  O(log  n)  time  using  n  processors  in  a  way 
which  is  described  after  Step  7.  For  this  it  is  crucial  to  observe  that  the  LCAs 
of  at  most  2n  edges  may  be  computed  as  part  of  this  Instruction. 

(d)  [PREORDER(uj^)  <  PREORDER(v)]AND[POSTORDER(uj^)  >  POSTORDER(v)].  Observe  that 
all  these  U},-s,  as  well  as  LCACv.u^  )  and  LCA(v,U£  ),  are  ancestors  of  v.  So,  they 
are  all  on  the  path  from  v  to  the  root  In  H  and  therefore  ANCESTOR(-,-)  is  a 
complete  order  on  them.  Find  the  lowest  ancestor  among  them  by  a  standard  minimum 
computation.  This  ancestor  is  the  required  Ad-LCA(v)  and  an  edge  that  caused  it 
is  the  required  e.  Complexitj^.  n  processors  are  used.  The  classification  into 
the  cases  takes  0(1)  time.  The  minimum  and  maximum  computations  in  (b),  (c)  and 
(d)  can  be  done  in  O(log  n)  time  using  the  obvious  algorithm  (that  is,  compare, 
Iteratively,  disjoint  pairs  of  "winners"). 

Step-4.-2.  Assign  to  each  edge  e  which  was  selected  at  Step  A.l  a  new  'serial 
number'  SERIAL(e)  as  follows.  SERIAL(e)  is  a  pair  which  consists  of  the  level  of 
Ad-LCA(v)  and  its  original  serial  number  (a  number  between  1  and  m).  That  is 
SERIAL(e)  :=  (LEVEL(Ad-LCA(V)) ."original  serial  number  of  e").  Complexity.  n 
processors,  0(1)  time.  Let  e  be  any  edge  in  E-T  such  that  SERIAL(e)  was  not 
defined  in  Step  4.2.  For  explanation  purposes  only  set 
SERIAL(e)  :=  (LEVEL(LCA(e)) ."origianl  serial  number  of  e"  +  m).  The  discussion 
below  explains  why  we  can  dispense  with  this  assignment  (and  actually  ignore  all 
edges  of  E-T  whose  SERIAL(e)  was  not  defined  in  Step  4.2  in  the  computation  of 
master  edges  that  follows)  without  affecting  the  correctness  of  the  algorithm. 

Let  e  and  f  be  two  edges  of  E-T.   Define  a   lexicographic  order   <^   between 
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SERIAL(e)  and  SERT.AL(f)  as  follows:  SERIAL(e)  <j  SERIAL(f)  if  SERIAL(ej).l  < 
SERIAL(E2).l  or  SERIAL(ej^ ) .  1  =  SERIAL(e2).l  and  SERIAL(ei).2  <  SERIAL(e2).2  , 
where  given  a  pair  (a,b),  (a,b).l  is  a  and  (a,b).2  is  b.  <•  is  a  complete  order  on 
the  SERIAL(-)  values  of  the  edges  of  E-T. 

Let  f  be  an  edge  of  T  .  Our  goal  Is  to  compute  MASTER(f).  It  is  the  edge  e  e E-T 
with  minimum  SERIAL(e)  whose  fundamental  cycle  contains  f. 

Claim.  Let  e=(u,v)  c  E-T  be  any  edge  whose  SERlAL(e)  was  not  defined  in  Step 
4.2.  It  is  impossible  that  e  is  the  master  edge  of  any  edge  of  T.  Proof.  Let 
w  -»^  u  be  the  father  edge  of  u  in  H.  If  (u,w)  is  included  in  the  fundamental  cycle 
of  (u,v)  then  due  to  the  fact  that  (u,v)  was  not  selected  in  Step  A. 2  there  must 
be  another  edge  which  is  adjacent  to  u  ,  say  (u,x),  such  that  LCA(u,x)  is  an 
ancestor  of  or  equals  to  LCA(u,v).  Therefore,  SERIAL((u,x))  <  SERIAL((u,v))  and 
(u,x)  "beats"  (u,v)  in  the  contest  to  become  the  master  edge  of  each  of  the  edges 
along  the  path  from  u  to  LCA(u,v)  in  H.  Similar  considerations  imply  that  (u,v) 
can  never  be  the  master  edge  of  any  edge  on  the  path  from  v  to  LCA(u,v)  in  H. 

Step  5.  Let  f=(u,v)  be  an  edge  of  T  and  say  that  u  is  the  father  of  v  in  H. 
For  each  such  edge  initialize  MASTER(f)  to  MIN{SERIAL(e)}  where  the  minimum  (with 
respect  to  <,)  is  taken  over  all  edges  e=(w,v)  in  E-T  whose  SERIAL(e)  is  defined 
in  Step  A. 2.  (There  are  <  n  such  edges).  Complexity.  n  processors,  O(log  n) 
time. 

The  interval  of  f  (eT)  in  the  Euler  path  of  T'  is  the  portion  of  the  Euler  path 
that  starts  at  the  u  -♦^  v  copy  of  f  and  ends  at  the  v  •+  u  copy.  Our  presentation 
identifies  MASTER(f )  with  both  MASTER(u  *  v)  and  MASTER(v  •♦  u).  No  confusion  will 
arise. 

Main  lemma.  Let  f=(u,v)  be  an  edge  in  T  where  u  is  the  father  of  v  in  H,  as 
before.  Then,  the  master  edge  of  f  is  the  edge  of  E-T  that  provides  the  minimum 
MASTER(-)  value  over  the  edges  in  the  interval  of  f.   Proof.    "Disconnect"  T  by 


"deleting"  f.  The  interval  of  f  in  the  Euler  path  of  T'  scans  the  component  of  v 
in  the  disconnected  tree.  The  master  edge  of  f  comes  from  the  edges  in  E-T  which 
have  one  endpoint  in  the  component  of  v  and  the  other  at  the  component  of  u.  All 
these  edges  are  taken  into  account  in  the  minimum  computation  for  f.  Let  e  be  an 
edge  of  E-T  whose  both  endpoints  are  in  the  component  of  v.  Observe  that  SERIAL(e) 
is  taken  into  account  in- the' minimum  computation  for  f.  This  seems  incorrect! 
However,  we  show  that  fortunately  SERIAL(e)  must  loose  in  the  minimum  computation 
to  an  edge  that  has  only  one  endpoint  In  the  component  of  f.  (This  explains  the 
role  of  our  refined  assignment  of  serial  number  to  edges  of  E-T).  There  must  be 
an  edge  g  in  E-T  that  has  only  one  endpoint  at  the  component  of  v  or  f  is  a  bridge 
of  the  input  graph.  LCA(g)  is  an  ancestor  of  v.  Thus,  LCA(g)  is  an  ancestor  of 
LCA(e),  and  therefore,  SERIAL(g)  <j  SERIAL(e). 

In  steps  6  and  7  we  compute  the  minima  described  in  the  main  lemma  for  every 
Interval  of  an  edge  In  T.  Assume,  w.l.g,,  that  the  length  of  the  Euler  path  (2n-2) 
is  a  power  of  two.  -    -^         -  ^  -   -   -•  -  ■ -^  ■■      ■-  - 

Step  6.  Compute  MIN{MASTER(-)}  over  all  edges  of  T'  (the  edges  of  the  Euler 
path).  For  reasons  that  will  become  clear  shorly  we  insist  on  computing  this 
mininmum  as  follows.  Apply  the  obvious  parallel  algorithm  which  uses  a  balanced 
binary  tree  whose  leaves  correspond  to  directed  edges  of  the  Euler  path  in  the 
order  of  the  path.  This  algorithm  finds  first  the  minimum  of  successive  pairs  of 
edges,  then  of  successive  fours,  then  of  successive  eights,  and  so  on.  This  can 
be  illustrated  by  an  upward  motion  in  Fig.  1:  Each  level  of  the  binary 
computation  tree  corresponds  to  a  stage  in  the  computation  of  Step  6.  Each  level 
of  Fig.  1  gives  the  segments  of  the  Euler  path  whose  minima  is  computed  in  the 
corresponding  stage  of  the  Step  6.  During  the  course  of  the  algorithm  we  save  all 
these  intermediate  minimum  computations.  Complexity.  n  processors,  O(log  n) 
time.   (Observe  that  n/log  n  processors  suffice). 
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Step  7.  The  interval  of  each  edge  f  of  T  can  be  represented  as  a  union  of 
O(log  n)  intermediate  minima  segments  which  are  pairwise  disjoint.  We  find  the 
minimum  of  SERIAL(-)  over  the  interval  of  f  by  computing  minimum  of  the  minima  of 
these  segments.  Complexity.  Assign  a  processor  to  each  of  the  n-1  edges  of  T.  It 
has  to  compute  which  segments  should  be  considered  and  compute  the  ralnlraura.  All 
this  needs  O(log  n)  time.  Next,  we  show  how  this  technique  can  be  used  to  compute 
the  lowest  common  ancestor  of  any  two  vertices  u  and  v  In  H.  (As  required  after 
case  (c)  of  Step  3).  Assume,  w.l.g.,  that  PREORDER(u)  <  PREORDER(v).  Consider 
the  Interval  of  the  Euler  path  that  starts  (exactly)  after  the  edge  of  the  form 
Wj  •♦  u  and  ends  before  the  edge  of  the  form  v  +  W2.  Observe  that  the  vertex  which 
yields  the  minimum  in  a  computation  of  minimum  level  over  the  vertices  of  this 
interval  is  LCA(u,v).  Levels  of  vertices  in  H  were  computed  earlier.  The  rest  of 
the  computation  of  LCA(u,v)  is  similar  to  steps  6  and  7. 

The  last  two  steps  of  the  algorithm  are  the  same  as  in  [At-82] . 

Step  8.  Assign  (arbitrarily)  a  direction  to  every  edge  in  E-T.  Complexity, 
m  -  (n-1)  processors,  0(1)  time.  For  the  alternative  implementation  observe  that 
this  can  be  done  by  p  processors  in  0(n^/p  +  1)  time.         -      :\^      _ 

Step  9.  To  every  (u,v)  e  T,  assign  a  direction  that  agrees  with  the 
fundamental  cycle  of  its  master.  This  is  done  as  follows:  Let  (i,j)  e  E-T  be  the 
master  of  (u,v)  and  assume,  w.l.g.,  that  u  is  the  father  of  v  in  H  and  that  (i,j) 
was  given  the  direction  i  -^  j  during  the  previous  step  of  the  algorithm.  One  of 
two  cases  can  occur:  (1)  If  v  is  an  ancestor  of  1,  then  the  u  +  v  direction  is 
assigned  to  (u,v).  (2)  If  v  is  an  ancestor  of  j,  then  the  v  +  u  is  assigned  to 
(u,v).   Complexity,   n  processors,  0(1)  time. 

Conclusion.   The  reader  is  invited  to  check  that  the  space  requirement  of  our  main 
inmplementation  is  0(n  +  m).   Most  of  our  use  of  the   doubling  procedure  was   in 


-m-  -  • 

order   to  rank  a  weighted   linked   list.    It   took  us  O(log  n)   time  using  n 

processors.   [Vi-8A]  presents  randomized  algorithms   that  run  almost-surely   in 

O(log  nlog*n)    (resp.    or   O(log  n))    time   using  n/(log  nlog  n)    (resp. 
nloglogn/log  n)  processors. 
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The  root  of  the  binary  computation  tree 
One  level  below  the  root 
Two  levels  below  the  root 
Three  levels  below  the  root 
Four  levels  below  the  root 
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